Relation between the High Density Phase and the Very-High Density Phase of 

Amorphous Solid Water 



^1" 
O 

o 

(N 



in 



Nicolas Giovambattista 1 , H. Eugene Stanley 1 and Francesco Sciortino 2 
1 Center for Polymer Studies and Department of Physics, 
Boston University, Boston, MA 02215 USA 
2 Dipartimento di Fisica and INFM Udr and Center for Statistical Mechanics and Complexity, 

Universtta' di Roma "La Sapienza" 
Piazzale Aldo Moro 2, 1-00185, Roma, Italy 
(Dated: 11 March 2004) 

It has been suggested that high-density amorphous (HDA) ice is a structurally arrested form of 
high-density liquid (HDL) water, while low-density amorphous (LDA) ice is a structurally arrested 
form of low-density liquid (LDL) water. Recent experiments and simulations have been interpreted 
to support the possibility of a second "distinct" high-density structural state, named very high- 
density amorphous (VHDA) ice, questioning the LDL-HDL hypothesis. We test this interpretation 
using extensive computer simulations, and find that VHDA is a more stable form of HDA and that 
in fact VHDA should be considered as the amorphous ice of the quenched HDL. 
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PACS numbers: 

The most common form of water in the universe is not 
a liquid but is a disordered solid named glassy water 
lj. Depending on the glass formation route — vapor de- 
position, hyperquenching of the liquid, pressure induced 
crystal amorphization — water forms amorphous solids of 
quite different structure, with the density varying by as 
much as 40 percent between the lowest density form and 
the highest 0, The structural and thermodynami- 
cal properties of the different amorphous forms of glassy 
water have been the focus of many recent experimental 
IlllJSJillJSIlSini, numerical and theoretical stud- 
ies [lj, UH' However, a clear picture of the phase 
diagram of glassy water is still missing. This is in part 
because the properties of glasses change drastically with 
aging|| 0| and because experiments at high pressure 
are difficult to perform. Glassy water can be found in at 
least two different forms, low-density amorphous (LDA) 
and high-density amorphous (HDA) ice |15l llq . If a 
glass is formed by extremely rapid cooling of a liquid, 
then one can naturally associate the glass with the liquid 
state. LDA can be obtained experimentally by very fast 
quenching of low-density liquid (LDL) water at normal 
pressure[l7j- However, the generation of an HDA glass 
via high pressure cooling of the liquid have proved an 
elusive goal; in fact, HDA is experimentally formed not 
by quenching but rather by compression either of LDA or 
of crystalline ice 0|. Therefore, the actual relationship 
of glassy water to the liquid at high pressure remains un- 
known. On the other side, computer simulations [TsUli^ 
indicate that HDA is a structurally arrested form of high- 
density liquid (HDL) water, and therefore, HDA can be 
obtained by quenching of HDL at high pressure. A sec- 
ond "distinct" high-density structural state, named very 
high-density amorphous solid (VHDA) has been recently 
discovered. The density of VHDA is 7-8% higher than the 
density of HDA. VHDA is generated by heating HDA un- 
der pressure. The resulting glass does not convert back 
to HDA when recovered at ambient pressure at T = 77K. 



The possibility of a VHDA phase in addition to an HDA 
phase raises many interesting questions po| . such as the 
relation between VHDA and HDA; this relation is im- 
portant to elucidate the hypothesis that below a critical 
temperature there are two distinct phases of liquid water, 
LDL and HDL 18]. Which phase, HDA or VHDA, more 
resembles liquid water at high pressure is particularly in- 
teresting, and it appears that the structural properties 
of VHDA are closer to liquid water at high pressure than 
those of HDA H H3- On the other side, the presence 
of two different high-density amorphous ices could imply 
the existence of more than one form of HDL as recently 
suggested by computer simulations plf. 

Computer simulations offer a potentially useful tool 
to probe the relation between HDA, VHDA and liquid 
water, since the time scale of most simulations is suffi- 
ciently short that liquid and glassy states can be studied 
for a wider range of state points than is accessible exper- 
imentally. Here we report a set of extensive simulations 
that suggest that VHDA not HDA, may be considered 
as a physical manifestation of the quenched high pres- 
sure liquid. Further, our simulations suggest that VHDA 
is not a new thermodynamically distinct structure but 
rather VHDA results from partial annealing of the HDA 
structures made possible by the higher annealing tem- 
perature. Hence HDA is not stable but rather is highly 
metastable, relaxing to VHDA in a fashion analogous to 
the way that, on slow heating, glasses generated with 
hyperquenched methods relax to glasses generated with 
standard cooling rates. 

Simulations offer the unique possibility of comparing 
(i) the glass resulting from conversion of HDA to VHDA 
with (ii) the glass generated by isobaric cooling of the 
liquid. However, a note of caution is in order, since in 
comparing experiments and simulations for glasses, one 
must carefully account for the significantly different time 
scales probed. State-of-the-art simulations probe time 
scales of 10-100 ns, so when the characteristic relaxation 
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time becomes comparable to this time, the system glas- 
sifies. Since the experimental homogeneous nucleation 
time is longer than 100 ns, glass configurations can be 
generated in simulations by cooling the liquid both at 
low and at high pressure. We shall exploit this fact in 
the present study. 

We simulate a system of 216 molecules, using the sim- 
ple point charge extended (SPC/E) model of water [2^ . 
This model has been studied extensively, and the p and 
T dependence of structural and dynamic properties in 
equilibrium are known. In particular, the SPC/E model 
reproduces the thermodynamic anomalies of water — e.g., 
it produces a maximum in p 23]. At p — 0.94 g/cm 3 at 
low T, the model describes well the LDA structure. We 
integrate the equation of motion using a time step of 1 fs 
and implement the reaction field method to account for 
long range forces. Results are averaged over 16 indepen- 
dent realizations. LDA ice configurations are generated 
by cooling to 77 K equilibrium liquid configurations at 
p = 0.90 g/cm 3 . To generate HDA configurations, we 
compress LDA at T = 77 K to P > 1 GPa at a rate of 
5 x 10~ 5 g/cm 3 /ps [3. The HDA to VHDA transition is 
studied heating HDA at constant pressure from T = 77 K 
using a heating rate of 30 K/ns. We have used the same 
heating/cooling and compression/decompression rates in 
all our calculations. 

FigureQ]compares, for a high pressure isobar, available 
experimental data (Fig. ^V) an d simulation results (Fig. 
^3). In Fig. the density p vs. temperature T for equi- 
librium liquid water is complemented with the p vs. T 
data of Ref. Q obtained by heating HDA from 77 K up to 
165 K, and then cooling back to 77 K. VHDA is the glass 
that results from densification of HDA during heating 
under high pressure. VHDA can be cycled, at constant 
pressure, between 77 K and 165 K without significant 
further density changes, leading us to hypothesize that 
VHDA is the more stable form of HDA. Our hypothe- 
sis is consistent with Mishima's observation that HDA 
samples annealed to 130-150 K at 1-1.5 GPa are charac- 
terized by identical scattering patterns [24|. Relaxation 
of HDA to VHDA, and the irreversibility of such transfor- 
mation upon further temperature-cycling, are facts rem- 
iniscent of slow heating of hyperquenched glasses, so the 
conversion of HDA to VHDA can be interpreted as a 
temperature-driven partial equilibration of the sample. 

In simulations HDA is generated — in analogy with the 
experimental procedure by compressing LDA or ice 1^ 
at 77 K [2{|, while LDA is generated by cooling of the 
liquid at ambient pressure. Figure ^3 shows p as HDA 
is heated. At about 100 K, p begins to increase, reach- 
ing at around 140 K a value 0.04 g/cm 3 larger than the 
density of HDA at 77 K. Decreasing T back to 77 K does 
not regenerate the HDA density, but instead p increases, 
in agreement with experimental results. The resulting 
denser material — which we identify with VHDA — can be 
cycled back and forth without significant changes in p. 



Next we test the hypothesis that VHDA is the glass 
that would be generated by cooling the high pressure liq- 
uid. To this end, Fig. [21 compares p for isobaric cooling 
of the liquid with p for VHDA, two systems generated by 
completely different thermal and pressure histories. The 
p data display remarkable similarity, supporting the pos- 
sibility that VHDA is indeed the physical realization of 
the glass generated by quenching the high pressure liq- 
uid. To further support this interpretation, Fig. |3| shows 
that the radial distribution functions of VHDA and of 
the glass obtained by isobaric cooling of the liquid are 
indeed indistinguishable. 

The results presented in Figs. |2 and indicate that 
compression of LDA at 77 K generates a system that 
is kinetically trapped due to the low temperature. This 
system relaxes to the more stable VHDA, whose structure 
is identical to the glass generated by quenching the high 
pressure supercooled liquid. In this respect, not HDA 
and LDA but rather VHDA and LDA should be thought 
of as the two distinct glassy states arising from the two 
distinct liquids associated with the hypothesized line of 
liquid- liquid phase transitions. 

If indeed HDA is a partially-equilibrated glass gener- 
ated by the compression technique, then one expects that 
VHDA should not convert to HDA under any transfor- 
mation which does not involve a different "intermediate" 
phase (such as LDA, or a crystal phase). To confirm this 
expectation, we recover VHDA at T = 77 K and ambient 
pressure with p = 1.26 g/cm 3 . Next we isochorically heat 
this system to 155 K (Fig.@}. We compare the structure 
of the resulting system with the structure of the VHDA 
glass recovered at the same density at T = 155 K. We 
find that the two glasses, which are generated from two 
completely different histories, are identical both thermo- 
dynamically (Fig. 4A) and structurally (Fig. 4B), and 
further we find that no VHDA to HDA transition takes 
place along the loop. Our conclusion motivates the need 
for a definitive experimental test that no VHDA — > HDA 
transition occurs. When VHDA is heated isochorically 
from 77 to 140 K, from a starting pressure of 0.02 GPa, 
the width of the x-ray pattern remains narrow like VHDA 
(and hence does not appear to revert to HDA); how- 
ever the position of its first maximum shifts (about 2 
degrees) in the direction of HDA (cf. Figs. 1A', IB', and 
ID' of Ref. 

In summary, our simulations suggest that VHDA is not 
a new form of glassy water but it is the result of annealing 
HDA upon heating at high pressure. We also find that 
VHDA is the glass obtained by fast quenching of high 
pressure water, i.e. HDL. How HDA and VHDA behave 
upon decompression to normal pressure and the phase 
diagram of amophous water are relevant topics that will 
be addressed in a separate work. 
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FIG. 1: Experimental data (from Ref. [jj) and simulation re- 
sults for the temperature dependence of the density p during 
the conversion of HDA to VHDA at high pressure (red line) 
and for the cycling of VHDA between 165 and 77 K (green 
and blue lines). The equilibrium p data in the liquid state 
(magenta line) from the HGK equation of state ]2& and for 
SPC/E potential are also shown. Numerical results for dif- 
ferent pressure values are qualitatively similar to the results 
shown in the figure for P=1.38 GPa. The density difference 
between HDA and VHDA decreases on increasing P. We also 
find that the density dependence of the VHDA recovered at 
T=77 K and ambient pressure varies between 1.22 and 1.28 
g/cm 3 . 
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FIG. 2: Comparison between the temperature-dependence of 
the density in VHDA (from Fig^B) and the temperature- 
dependence of the density during continuous cooling of the 
liquid down to 77 K, at a quenching rate of — 10 4 K/ns. 
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FIG. 3: Comparison between the structural properties of 
VHDA (green) and the glass obtained by cooling the liq- 
uid under pressure calculated at T = 77 K and P — 1.38 
GPa (magenta). Both oxygen-oxygen, oxygen-hydrogen, 
hydrogen-hydrogen radial distribution functions are shown. 
Despite the extremely different previous histories, the VHDA 
structure appears identical to the structure of the liquid 
cooled under pressure. 
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FIG. 4: Demonstration that at T = 155K and p = 1.26 g/cm 3 
two glasses, each with completely different histories (red and 
blue curves), are identical both (A) thermodynamically and 
(B) structurally. (A) P — p diagram. Blue curve: recovering 
of the VHDA at ambient pressure at T = 77K followed by 
an isochoric heating at density 1.26 g/cm up to T — 155-K". 
Red curve: isothermal decompression of VHDA at T = 155K 
down to density 1.26 g/cm 3 . The final product is indepen- 
dent of the path despite the different histories. The impor- 
tant point is that no transformation from VHDA to HDA is 
observed. (B) Comparison between the radial distribution 
function of the two glasses obtained following the blue and 
the red paths in (A). 



